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ABSTRACT 



o 

1-^ ■ Active Galactic Nuclei (AGN) consist of an accretion disk around a super- 

I ' massive black hole which is in turn surrounded by an obscuring torus of dust 

;_i ■ and gas. As the resulting geometry of this system affects the observable proper- 

^ ■ ties, quantifying isotropic indicators of intrinsic AGN luminosity is important in 

selecting unbiased samples of AGN. In this paper we consider five such proxies: 

the luminosities of the [OIII]A5007 line, the [OIV]25.89/xm line, the mid-infrared 

(MIR) continuum emission by the torus, and the radio and hard X-ray (E > 



O . lOkeV) continuum emission. We compare these different proxies using two com- 

Q ■ plete samples of low-redshift type 2 AGN selected in a homogeneous way based 

f-~-. ■ on different indicators: an optically selected [OIII] sample and a mid-infrared 



selected 12/im sample. To assess the relative merits of these proxies, we have 
undertaken two analyses. First, we examine the correlations between all five 
different proxies, and find better agreement for the [OIV], MIR, and [OIII] lumi- 
S^ . nosities than for the hard X-ray and radio luminosities. Next, we compare the 

ratios of the fluxes of the different proxies to their values in unobscured Type 
1 AGN. The agreement is best for the ratio of the [OIV] and MIR fluxes, while 
the ratios of the hard X-ray to [OIII], [OIV], and MIR fluxes are systematically 
low by about an order-of- magnitude in the Type 2 AGN, indicating that hard 
X-ray selected samples do not represent the full Type 2 AGN population. In a 
similar spirit, we compare different optical and MIR diagnostics of the relative 
energetic contributions of AGN and star formation processes in our samples of 
Type 2 AGN. We find good agreement between the various diagnostic parame- 
ters, such as the equivalent width of the MIR polycyclic aromatic hydrocarbon 
features, the ratio of the MIR [OIV] /[Nell] emission- lines, the spectral index of 
the MIR continuum, and the commonly used optical emission-line ratios. Finally, 
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we test whether the presence of cold gas associated with star-formation leads to 
an enhanced conversion efficiency of AGN ionizing radiation into [OIII] or [OIV] 
emission. We find that no compelling evidence exists for this scenario for the 
luminosities represented in this sample (L^o/ ^ 10^ - 8 x 10^^ Lq). 



1. Introduction 

Active galactic nuclei (AGN) are powered by accretion onto a central supermassive 
blackhole. According to the unified model (e.g. Antonucci 1993), an optically-thick torus of 
dust and gas surrounds this central engine. Hence orientation of the system plays a central 
role in determining the observable features of AGN. In Type 1 AGN, the system is oriented 
face-on, leaving an unobstructed view of the central engine and the broad line region. In 
contrast. Type 2 AGN are oriented edge-on, blocking the accretion disk and the broad line 
region. These obscured AGN can be identified by their narrow optical and mid-infrared 
emission lines which originate in gas photoionized by accretion disk photons. This narrow 
line region (NLR) extends hundreds of parsecs away from the central source and is therefore 
not significantly affected by torus obscuration. 

The luminosity of emission lines formed in the narrow line region can therefore be 
used as isotropic indicators of intrinsic AGN luminosity. The fiux of the [OIII]A5007 line 
is commonly used as such a diagnostic (e.g. Bassani et al. 1999, Heckman et al. 2005) as 
it is one of the most prominent lines and suffers little contamination from star formation 
processes in the host galaxy. This line can be attenuated by dust in the host galaxy, though 
this effect can be somewhat remedied by applying a reddening correction using the observed 
Balmer decrement (i.e. the observed ratio of the narrow Ha/ H/3 emission-lines compared t o 



the intrinsic ratio) and the extinction curve for galactic dust (jOsterbrock fc Ferland II2006I ). 



Isotropic indicators of AGN luminosity also exist in the infrared band and are much less 
affected by dust extinction than the optical [OIII] line. Recently, the luminosity of the [OIV] 
25.89/im line has been shown to be a robust proxy of AGN power (e.g. Melendez et al. 2008, 
Rigby et al. 2009, Diamond-Stanic et al. 2009): it is formed in the NLR, so it is not affected 
by torus obscuration, and with an ionization potential of 54.9 eV, starburst activity does 
not significantly contribute to this line. AGN also emit over 20% of their bolometric fiux 
in the mid- infrared (MIR), where photons produced by the continuum are absorbed by the 
torus and re-radiated (e.g. Spinoglio & Malkan 1989). This MIR emission from the dusty 
torus can also be a proxy for the intrinsic AGN luminosity. Two potential issues with the 
MIR continuum are contamination by emission from dust heated by stars (e.g. Buchanan et 
al. 2006; Deo et al. 2009) and possible anisotropy in the torus emission (e.g. Pier & Krolik 
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1992, Buchanan et al. 2006, Elitzur & Shlosman 2006, Nenkova et al. 2008). 

Radio and hard X-ray (E > 10 keV) flux can serve as proxies of the intrinsic AGN con- 
tinuum. Radio emission has been shown to be similar between type 1 and type 2 AGN (e.g. 
Giuricin et al. 1990, Diamond-Stanic et al. 2009, Melendez et al. 2 010) and correlated with 



optical luminosity, in particularly the [OIII] flux (jXu et al. Ill999l ). making high resolution 
radio observations that isolate emission from the nucleus another diagnostic of intrinsic AGN 
power. Hard X-rays can pierce through the obscuring torus, provided that the object is not 
heavily Compton thick (Nh < 10^^ cm~^), and has therefore been used as a method to select 
AGN samples (e.g. Winter et al. 2008, Treister et al. 2009). 

Connections have been observed between star formation activity in the host galaxy and 
the central AGN (e.g. Kauffmann et al. 2003). Starburst activity can be parametrized by 
various IR features, such as the equivalent width (EW) of polycyclic aromatic hydrocarbons 
(PAHs). IR and optical data, such as the ratio of fine structure lines and the shape of the 
spectral slope, can also reveal the relative amount of AGN to starburst activity. Samples of 
Seyfert 2 galaxies (the predominant local class of type 2 AGN) are useful in examining the 
relationships among these star formation vs. AGN-to-starburst indicators as the obscuration 
of the central engine allows detailed study of the host galaxy. 

To address these issues, we will use two complete and homogeneous Sy2 samples selected 
based on isotropic indicators of AGN luminosity (one [OIII]-selected sample and one MIR- 
selected sample). We will compare the various diagnostics of intrinsic AGN luminosity and 
probe for biases resulting from sample selection criteria, starburst contamination, errors 
introduced from extinction correction, and scatter due to the various physical mechanisms 
producing these emission features. Such biases are likely minimized in the diagnostic ratios 
with the smallest dispersion. Where available, we compare these ratios with the Syl values to 
probe for differences due to the inclination of the system, thus testing to which extent these 
indicators of intrinsic AGN luminosity are truly "isotropic." We will also test the agreement 
among mid-infrared and optical star formation indicators. Finally, we will examine the 
possibility that the fraction of the AGN ionizing luminosity that is converted into [0111] and 
[OIV] emission is systematically higher in systems in which there is a copious supply of dense 
gas associated with starburst activity. 
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2. The Data 

2.1. Sample Selection 

The selection of the SDSS [OIII] sample is discussed in detail in LaMassa et al. (2009). 
In brief, Type 2 AGN were drawn from a parent sample of approximately 480,000 galaxies in 
SDSS Data Release 4 = DR4 (Adelman-McCarthy et al. 2006). The Type 2 AGN with an 
observed [OIII] flux greater than 4 x 10~^^ erg s~^ cm~^ were selected, providing a complete 
sample of 20 Sy2s (hereafter the "[OIII] sample," listed in Tabled]). 

The mid-IR sample com prises the Seyfert 2 galaxies from the original IRAS 12yum survey 



( jSpinoglio &: Malkan I Il989l ). This represents a complete sample of Sy2s down to a fiux- 
density limit of 0.3 Jy at 12 /im, drawn from the IRAS Point Source Catalog (Version 2), 
with latitude \b\ > 25° to avoid contamination from the Milky Way. We have dropped NGC 
1097 from this original sample as it has since been classified as a Type 1 Seyfert (Storchi- 
Bergmann et al. 1993), leaving 31 mid-IR selected Sy2s (hereafter the "12/im sample," listed 
in Table |2]). 



2.2. Optical Data 

The optical data for the [OIII]-sample were drawn from SDSS DR4, whereas the optical 
data for the 12/im sample were collected from the literature or from SDSS Data Release 
7 (DR7) where available. The reddening corrected [OIII] flux (Fi^oiii],corr) was calculated 
using the observe d Ha /H/3 ratio and an i ntrin sic ratio of 3.1 with R=3.1 extinction curve 



for galactic dust f lOsterbrock fc Ferland 1 120061 ). Tables [3] and H] list the optical emission 



line fluxes and ratios utilized for this study, as well as the relevant literature sources for 
the 12/im sample. The black hole masses (Mbh) were derived for the [OIII] sample by the 
SDSS velocity dispersion (a) and the M-a relation {Mbh = 10^-^3(cr/200 km s"^)^-"^ Mq, 
Tremaine et al. 2002). We used literature values for Mbh for the 12/im sources, with most of 
the masses derived using the M-a relation cited above. For F04385-0828, F05189-2524 and 
TOLOLO 1238-364, the full width half max (FWHM) of the [OIII] line was used as a proxy 
for the velocity dispersion (Wang & Ahang 2007, Greene & Ho, 2005), and photometry of 
the host galaxy was used to estimate Mbh for F08572+3915 (see Veilleux et al. 2009 for 
details) . 
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2.3. Infrared Data 

The infrared data presented here were obtained from the Infrared Spectrograph [IRS, 
Houck et ah 2004) on board the Spitzer Space Telescope. Low-resolution spectra were ob- 
tained using the Short-Low (SL, 3.6" x57" aperture size) and the Long-Low (LL, 10.5" xl68" 
aperture size) modules and high-resolution spectra were provided by the Short-High (SH, 
4.7"xll.3" aperture size) and the Long-High (LH, lLr'x22.3") modules. 

The Sy2s in the [OIII]-sample were observed in IRS staring mode in both high and low 
resolution under Program ID 30773. For the 12-/im sample, high resolution data existed for 
all 31 Sy2s but low resolution data were only available for 30 galaxies (IRAS 000198-7926 
lacked low resolution data). The high resolution data were obtained in IRS staring mode 
for 30 of the Sy2s (NGC 5194 had only IRS spectral mapping mode high-resolution data in 
the archive). Several galaxies had multiple IRS observations: we analyzed these observations 
independently and compared our results between the two observations. For the low-resolution 
data, IRS staring mode was used when available with the remainder observed in IRS spectral 
mapping mode. The Spectral Modeling Analysis and Reduction Tool (SMART, Higdon et 
al. 2004) was used to reduce the staring mode observations for the 12yU,m-sample to be 
consistent with previous IRS analysis of the 12 fun sample (e.g. Tommasin et al. 2008, Wu 
et al. 2009, Buchanan et al. 2006), Spitzer IRS Custom Extraction (SPICE) was used to 
analyze the staring mode observations for the [OlllJ-samplqj, and the Cube Builder for IRS 
Spectra Maps (CUBISM, Smith et al. 2007a) was utilized to analyze the spectral mapping 
observations. Table [5] lists the Program ID(s) for each galaxy, the IRS mode used, and 
spectral extraction area for low- resolution spectral mapping mode data (discussed below). 



2.3.1. High- Resolution IRS Staring Spectra 

We used the basic calibrated data (BCD) pipeline products as the starting point for 
our analysis. Rogue pixels were removed using the IDL routine IRSCLEAN_MASK and the 
rogue pixel mask matching the campaign number of the observation. Dedicated off-source 
background observations were taken for all sources in the [OIII]-sample and for most of 
the Sy2s in the 12/im sample. Multiple background observations, if present, were coadded 
within each nod and subsequently subtracted from the source image. The background- 
subtracted source images were then coadded between the two nods. The galaxies in the 12/xm 



^Though the IRS staring data for the 12/im sample was reduced using SMART and the [OIII]-sample with 
SPICE, the effect of the different reduction software is expected to be ncghgible on the derived parameters 
used in this analysis. 
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sample that had dedicated background observations and thus were background subtracted 
are marked with a "b" in Table [51 For sources in the 12/ini sample without dedicated 
off-source observations, no background subtraction was performed. 

Spectra were extracted from these combined observations, using the full aperture ex- 
traction mode. The edges of each order were then inspected, removing any data points that 
fell outside of the calibrated range for that order (IRS Data Handbook, Version 3.1, Table 
5.1). The orders were then combined using a 2.5-cr clipping mean, resulting in a final cleaned 
spectrum. 



2.3.2. Low- Resolution IRS Staring Spectra 

The low-resolution data were processed in a similar fashion as the high-resolution data, 
i.e. we started with the BCD products and removed the rogue pixels with IRSCLEAN_MASK. 
However, for these observations a background data set was built for each nod and order by 
coadding the off-source order and nod position. The background-subtracted nods (following 
the same procedure as above) were combined for each order and the spectra then extracted 
using tapered column extraction. The orders were combined using a 2.5-a clipping average. 
This procedure was executed separately for the SL and LL module. Fourteen of our galaxies 
had low resolution IRS staring mode data; the rest were acquired in spectral mapping mode. 
We note that IRAS 00198-7926 did not have archival low resolution spectral data. 



2.3.3. Spectral Mapping Spectra 



The IRS spectral mapping observations were analyzed with CUBISM (jSmith et al. 



2007al ). which uses the BCD data to create 3-D spectral cubes (one spectral dimension and 
2 spatial dimensions). For the low resolution data, background observations were built from 
the other order of the on-source module (e.g. SL 2 was used as the background for SL 1, 
etc.). After the rogue pixels were removed, using the default "autogen bad pixels" option in 
CUBISM, a spectral cube was built. Spectra were then extracted using matched apertures 
among the detectors and centered on the nucleus. The aperture extraction size for these low- 
resolution spectral mapping observations are listed in Table [51 The low resolution spectral 
mapping data for NGC 1068 was saturated near the nucleus and consequently not included 
in this analysis. 

For the IRS spectral mapping high resolution observation of NGC 5194, no background 
subtraction was performed. The spectrum was extracted over the full cube, corresponding 
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to a size of 31.5" x45" in the LH module and 13.8" x 27.6" in the SH module. 

2.4. Radio and Hard X-ray Data 

The radio and hard X-ray data were drawn from the literature; VLA radio data at 
8.4 GHz were only available for the 12/im sample (Thean et al. 2000). In several cases, 
multiple radio components were analyzed; we included only the flux for the component that 
was nearest the published center of the galaxy. Twenty-six of the 31 12/im sources had radio 
data, with 3 additional sources having upper limits. The ha rd X-ray fluxes origina ted from 



the 2 2-month Swift-BAT Sky Survey (JTueller et al. 1 12010| ) and from BeppoSax (JDadina 



20071 ). Only 11 out of the 31 12//m sources and one of the 20 [OIII] sources (IC 0486) have 
X-ray detections in the 14-195 keV range. We adopted an upper limit of 3.1x 10~^^ erg cm~^ 
s~^, the flux limit of BAT, for the remainder of the sample when an upper limit was not 
quoted in either Tueller et al. (2009) or Dadina et al. (2007). 



3. Measurements 

3.1. IR Emission Line Fluxes 

The high resolution spectra were utilized to measure the emission line fluxes: a Gaussian 
profile was fit to the emission line feature, with the local continuum, centered on the line's 
rest-frame wavelength, fit by a zero- or first-order polynomial. The errors were estimated 
by calculating the root-mean-square (RMS) around this local continuum and measuring the 
fiux values with the continuum shifted by ± the RMS. In the cases where an emission line 
was not present, a 3-cr upper limit was estimated from the RMS around the best-fit local 
continuum (where the RMS is assumed to be the 1-a error). In the cases with multiple 
observations per galaxy, we measured the emission line fiuxes independently and averaged 
the resulting values; these fiux measurements agreed within several percent between most 
of the individual observations, with at most a factor of ~1.5 discrepancy, which was only 
present in one of the sourceso Tables [6] and [7] list the emission line fiux values for the 
[OIII] and 12/im samples, respectively. Comparing our line fiux values with Tommasin et 



^We note that NGC 1143/4 has two high rcsohition archival observations, one centered at RA=43.8004, 
Dec=-0.1839, and the other at RA=43.7985 and Dec=-0.1807, a distance of -13.4". We present the line 
fluxes from the first region as this corresponds to NGC 1144 which is classified as a Sy2 in SIMBAD. The 
optical data are for NGC 1144 as well. 



al. (2008, 2010), we find that our [OIV] flux values largely agree within a factor of 1.5 
(with the exception of NGC 1667 and NGC 7582 where their values are a greater than a 
factor of 2 higher than ours). However, their [Nell] flux values are generally systematically 
higher by a factor of ~2.5 - ~4.5, though we do obtain consistent values for NGC 424, 
NGC 5135 and NGC 5506. Despite these differences in the measured [Nell] line strength, 
we obtain similar results to Tommasin et al (2010), namely that as the relative contribution 
of the AGN to the ionization field increases (parameterized by [OIV] /[Nell]), the starburst 
strength (parameterized by the PAH equivalent width) decreases. 



3.2. IR Continuum Flux and PAHs 

The MIR continuum flux values (Fmir) and PAH equivalent widths (EWs) were mea- 
sured using the low resolution spectra. For the galaxies that had multiple observations, we 
utilized the observations that had consistent flux values in the overlap region between the 
SL and LL modules: Program ID 30572 for For NGC 1386, NGC 4388, NGC 5506 and NGC 
7130; Program ID 0086 for NGC 5135; and Program IDs 00086 and 30572 for NGC 5347 (for 
this source, the analysis was done separately for each observation and the results averaged 
together). The MIR continuum flux was measured at 13.5 yum (rest-frame), averaged over a 
3/im window; these flux values are listed in Table E] for the [OIII]-selected sample and Table 
[7] for the 12/im sample. This window was chosen as it is free from strong emission line and 
PAH features jj In LaMassa et al. (2009), we included the flux centered at 30/xm as part 
of the MIR flux diagnostic. However, emission at this higher wavelength can be strongly 
affected by star formation processes in the host galaxy (e.g. Deo et al. 2009, Baum et al. 
2010), so here we use Fias^^ as Fmir- 



We used PAHFIT (ISmith et al Il2007bl ) to measure the PAH EWs, a program which uses 
a model consisting of several components: a starlight component represented by blackbody 
emission at T = 5000 K, a thermal dust continuum constrained to default temperature bins 
(35, 40, 50, 65, 135, 200 and 300 K), IR emission lines, PAH (dust) features and extinction 
(we used a foreground extinction screen). As PAHFIT requires a single stitched spectrum, 
the SL spectrum was scaled to match the LL spectrum, with typical adjustments under 20% 
(though several galaxies were adjusted by ~40% and NGC 7582 by greater than a factor of 
6, indicating the presence of extended IR emission in this object). Here we utilize the EW of 
the PAH features at 11.3//m and 17/im, which consist of the features within the wavelength 



■^Though this range docs include the [Nell] 12.81/im Hnc, in most cases this comprises less than 1% of 
the MIR flux, with the exception of NGC 7582 where the [Nell] line is ~1.5% of the MIR flux. 
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range 11.2-11.4yuni and 16.4-17.9/im, respectively (Tables |6] and [7]). However, we note that 
the current version of PAHFIT has a bug which assumes the PAH EW feature to be Gaussian 
rather than Drude, which could underestimate the PAH EW by a factor of 1.4. We report 
the EWs as reported from PAHFIT, with the caveat that these may be lower limits. 

We compared our results with the 11.2/im feature from Wu et al. (2009) and the 
11.3/im and 17yU.m features measured by Gallimore et al. (2010), where in the latter, we 
added their published 11.2yum and 11.3/im EW values. Wu et al. employed a spline fit 
between 10.8/im and 11.8/im to measure the EW. With this method, the results are widely 
influenced by the choice of anchor points for fitting the pseudo-continuum and can result 
in an underestimat e of the EW compa red to a method that utilizes spectral decomposition. 



such as PAHFIT fISmith et al I l2007bl ). Of the 28 sources we have in common with Wu 
et al. (2009), 12 of them had consistent 11.3/im EW values (within a factor of 2), 6 had 
lower values than we obtained (which would be expected from the disagreements between 
the spline vs. decomposition methods mentioned above) and 10 had higher values, where 
for 6 of these, PAHFIT had obtained an EW value of zero, yet the spline method yielded 
a measurement. Comparing our results with Gallimore et al. (2010) gave better results, 
though a discrepancy did still exist: of the 23 sources in common, we obtained consistent 
EW values (within a factor of 2) for 11 sources at 11.3/im and for 12 sources at 17/im. Though 
Gallimore used PAHFIT to measure these features, they modified the code to include more 
fine-structure lines, fit silicate emission features, and use the cold dust model from Ossenkopf 
et al. (1992); they also generated their own software to build spectral data cubes whereas 
we employed CUBISM. Such differences could account for the inconsistencies in our PAH 
EW measurements. Though the derived EWs are different from those reported by Wu et al. 
and Gallimore et al. for at least half the sources we have in common, our main conclusions 
based on PAH EWs agree qualitatively with Wu et al. and Gallimore et al.: PAH features 
are associated with other star formation activity indicators (Gallimore et al. 2010, Wu et 
al. 2009) and the EWs are inversely correlated to the strength of the ionization field (Wu et 
al. 2009 where they use the IRAS colors to parameterize AGN strength). 

In the discussion that follows, we divide the 12/im-sample into two classes, those with 
weak PAH emission ("PAH- weak" sources) and those with strong PAH emission ("PAH- 
strong" sources, galaxies with EW > 1 /im in either the 11.3 /im or 17/im band, with PAH 
EWs detected in both bands); the strong PAH emission is likely due to starburst activity in 
the host galaxy (see §5.2). 
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4. Diagnostics of Intrinsic AGN Luminosity 

Our goal is to evaluate the relative efficacy of the five different proxies for the AGN 
intrinsic luminosity under consideration in this paper. We expect that these different proxies 
will not agree perfectly, due to the different physical mechanisms that produce and affect the 
emission features as well as biases resulting from sample selection, starburst contamination, 
statistical errors and in some cases, uncertain application of extinction corrections. To 
address this, we will undertake two kinds of comparions. 

First, we will use our two Sy2 samples to inter-compare these proxies in a pair- wise 
fashion and measure the amount of scatter in the corresponding flux ratios. Which proxies 
agree best with one another? Second, we will compare these pairs of flux ratios to the 
corresponding values for unobscured Type 1 AGN to test which proxies are more "isotropic," 
i.e. suffer the least AGN-viewing-angle dependence. 

Figures ID - [TUJ show the histograms of a subset of ratios for the flve proxies. In each plot, 
the solid black line represents both samples combined, the red dashed line and green dotted- 
dashed line delineate the 12/im sample ("PAH-weak" and "PAH-strong" sources respectively) 
and the cyan flUed histogram reflects the [OIII]-sample. Adjacent to these histograms are 
the luminosity vs. luminosity plots, showing the correlation between these indicators: the 
cyan asterisks represent the [OIII] sample, the red diamonds (green triangles) depict the 
"PAH-weak" ("PAH-strong") 12/im sources, and the dashed black line represents the best 
flt from multiple linear regression analysis (i.e. the REGRESS routine in IDL), where in the 
flgure captions, p is the linear regression coefficient and Puncorr is the probability that the two 
quantities are uncorrelated. Though the distance dependence in luminosity vs. luminosity 
plots enhances the correlation compared to flux vs. flux plots, we employed this method as 
the 12/im sample lies at a systematically lower redshift, and thus have higher flux values, 
than the [OIII] sample. One of our main goals is to examine the dispersion in the flux ratios, 
where this distance dependence cancels out. In §4.3, we test if these ratios are affected by 
luminosity. 

Where available, the values for Syls are included in these plots. The results are sum- 
marized in Tables [8] and [9] which lists the mean and sigma of each ratio for the combined 
sample and the sub-samples separately. In the histograms and luminosity plots, the upper 
limits are plotted but not included in the analysis of the mean and sigma (except for the 
ratios involving the hard X-ray flux). Since only 12 of the 51 AGN were detected in hard 
X-rays, we have employed survival analysis to quantify the correlations among the proxies 
and to calculate the mean of the ratios. This approach takes the upper-limits into account 
(ASURV Rev 1.2, Isobe and Feigelson 1990; LaValley, Isobe and Feigelson 1992; for uni- 
variate problems using the Kaplan-Meier estimator, Feigelson and Nelson 1985; for bivariate 
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problems, Isobe et al. 1986). 



4.1. Inter-Comparison of Proxies 

The isotropic luminosity diagnostics tliat agree best, and therefore may be least subject 
to the uncertainties and errors discussed above, are F[o///],ofes5 F[o/y] and Fmib- A wider 
spread is present between the radio and hard X-ray fluxes compared with the optical and 
MIR values. 

In all cases, a wider dispersion is present between all the flux ratios in the 12/im sample 
as compared to the [OIII] sample. Below, we examine whether such a scatter could be due 
to aperture effects, extinction corrections applied to the [OIII] flux, starburst contamination 
to the the MIR flux or if it represents a real difference between AGN selected on the basis 
of [OIII] flux versus MIR flux. 

Since the 12/im Sy2s are typically more nearby than the [OIII]-selected galaxies, aperture 
effects can potentially play a signiflcant role when comparing flux values by either missing 
NLR flux or obtaining too much host galaxy contamination. However, we flnd no evidence 
in our data for such an effect (see Appendix). Another possible explanation for this wider 
dispersion is that the optical data for the 12/im sample are drawn from the literature, which 
can introduce scatter into the optical diagnostics as such data are not taken and reduced 
in a uniform matter. The most striking example of this is the comparison of the [OIV] flux 
with the observed and extinction corrected [OIII]-flux: de-reddening the [OIII]-flux widens 
the dispersion in the 12/im-sample (see Figured]). This can be an artifact of using literature 
Ha and H/3 values and to a lesser extent can be due to large amounts of dust in the host 
galaxies of the 12/im sample as evidenced by the wide range of Balmer decrements. Gould- 
ing and Alexander (2009) note that galaxies that would not be identifled optically as AGN 
(i.e. have a low "D" value, see §5.1) tend to have similar Balmer decrements yet higher 
F[o/v]/F[o///],o6s ratios than optically identifled AGN. This result suggests that applying a 
reddening correction using the Balmer decrement still under-represents the intrinsic [OIII] 
flux. However, our 12/im sources with lower "D" values (<1.2) do not show systematically 
higher F[oiv]/Fioiii],obs ratios (see Figure [11]), indicating that "extra" extinction that Gould- 
ing and Alexander observe in their sources is not present in ours. As expected, the ratio of 
F[oiv]/F[oiii],obs increases with Ha/H/3, denoting that both quantities trace host galaxy ex- 
tinction, though with wide scatter. Comparison with the locus of points for the [OIII] sample 
shows that the Balmer decrement is systematically higher for similar 12/im Fioiv]/F[oiii],obs 
values, suggesting that the the 12/im Balmer decrements are over-estimating the amount of 
dust present rather than under-estimating. This result indicates that the literature Balmer 
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decrements, which are not analyzed in a systematic and homogenous way, are introduc- 
ing uncertainties that bias the results and do not better recover the truly intrinsic AGN 
luminosity. 

However, this can not be the only cause of the greater scatter in the 12/im sample, since 
the ratio of MIR/[OIV] fluxes shows more scatter in the 12 fun sample (Figure [3]), though 
these data were analyzed homogeneously. Could the presence of Sy2s in the 12yum sample 
that have significant contributions from starburst activity create the wider dispersion in these 
diagnostics? To address this issue, we isolated the "PAH-strong" sources, which have greater 
amounts of star formation activity (discussed in detail below). The distributions between the 
"PAH-strong" and "PAH-weak" sub-samples are similar, suggesting that starburst processes 
are not responsible for the wider dispersion. We also focused on sources with a limited 
"D" value, which as noted above indicates the relative contribution of AGN to starburst 
activity. Repeating the calculation of mean and standard deviations on the flux ratios for 
the 12/im sources with 1.2< D < 1.7 did not result in a significant decrease (a factor of 
2 or more) in the dispersion with the exception of log (FMiR/F[oiii]^obs) ((7=0.42 dex), log 
(Ppiiiiobs/FsAGHz) (o-=0.60 dex) and log {F[oin],corr/FsAGHz) (cr=0.51 dex). For these first 
two ratios, this is largely due to the removal of the 3 outliers (F04385-0828, F08572+3914 
and Arp 220) with systematically higher (lower) Fm/r/F[o///],o6s {F[oiii],obs/F8AGHz) values 
from the full sample. The dispersions for the other ratios were still systematically higher 
than the [OIII]-sample. We conclude that there is a real difference between the AGN selected 
on the basis of [OIII] emission-lines and MIR continuum. 

We also compared our results with two other samples of Seyfert 2 galaxies: one a 
complete sample down to a flux limit of (1-3) xlO^^^ erg cm"^ s~^ at 14 - 195 keV drawn 
from the 3 and 9 month Swift-BAT survey (Melendez et al. 2008 and references therein) and 
the other drawn from the revised Shapley-Ames catalo g (Shapley fc Ames 1932; San dage & 
Tammann 1987), consisting of galaxies with Bt < 13 ( Diamond-Stanic et al. II2009I ). Here 



we include those radio quiet Seyfert types 1.8-2 that have measured [OIII] and [OIV] fluxes, 
giving 12 and 56 Sy2s, respectively. The log of the ratios of [OIV] to [OIII],obs for both 
samples are higher than our combined sample, 0.60 ± 0.74 dex (Melendez et al. 2008) and 
0.57 ± 0.67 dex (Diamond-Stanic et al. 2009) vs. 0.08 ± 0.41 dex, but the differences are 
not statistically significant. A wider dispersion is present in these comparison samples than 
the [OIII]-selected sample (as was the case for the 12/im sample). This may indicate that 
selection based on [OIII] leads to better agreement between between the [OIII] and [OIV] flux 
rather than selection based on other methods. This effect could be due to the [OIII]-bright 
sources having less extinction in the NLR than Sys selected in other ways. 

As the samples in Weaver et al. (2010) and Winter et al. (2010) samples were selected 
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based on their hard (14-195 keV) X-ray flux from the Swift-BAT 9 month (Winter et al. 
(2010) and Weaver et al.(2010)) and 22 month (Weaver et al. (2010)) catalog, we can 
compare our Sy2 X-ray flux ratios. Using the values in Winter et al. (2010), we flnd log 
{Fu-i95kev/'P[oiii],obs) = 2.76±0.59 dex and log {Fi4_iQ5kev/F[oiii],corrP = 2.34± 0.69 dex 
for Sy2 galaxies. The log (Fi4_i95fce\//F[o/y]) ratio from Weaver et al. (2010) for Sy2s 
is 2.38±0.45 dex. All three values are systematically higher than what we obtain for our 
samples of Sy2 galaxies by roughly an order of magnitude (see Table [9]). This is depicted 
graphically in Figure [T21 Employing survival analysis, we compared these ratios between 
the BAT-selected Sy2s and the [OIII] and 12//m samples separately and flnd that they differ 
signiflcantly (i.e. p <0.05, corresponding to the 2a level, that they are drawn from the 
same parent population), with the caveat that with only one [OIII] Sy2 detected by BAT, 
the comparison between BAT and [OIII] selected Sy2s may be less robust. These differences 
suggest that the samples selected in hard X-rays do not fairly sample the population of Type 
2 AGN selected in the MIR and possibly the optical, however comparisons with BAT-selected 
Syls reveal mixed results (see §4.2). 



4.2. Comparison with Syls 

In order to determine if the proxies we are considering are affected by the orientation 
of the AGN, and evaluate the extent to which they may be considered truly "isotropic," we 
compared our results for Sy2 with the corresponding values for Syl, using data taken from 
the literature. The Syl MIR fluxes were calculated from the 14.7/xm flux densitites reported 
in Deo et al. (2009), where they analyzed a heterogeneous sample of Syl and Sy2 galaxies 
available in the Spitzer archive, ranging in redshift 0.002 < z < 0.82. The radio flux values 
were derived from the high-resolution 8.4-GHz flux density values from Thean et al.(2000), 
which presented analysis of the extended 12/im sampleo The hard X-ray data (14-195 keV) 
are drawn from Mel endez et al. (2008, sample selection described above), the 22-month 



Swift-BAT Catalog (jTueller et al. Il2010[ ) and Rigby et al. (2009, same parent sample as 
Diamond-Stanic et al. 2009, with X-ray data derived from the 22-month Swift-BAT Catlog, 
BeppoSAX (Dadina 2007) and Integral (Krivonos et al. 2007)). The comparison Syl [OIII] 



"'We note that for the cases where the Winter et al. 2010 Balmer decrement was less than the assumed 
intrinsic value (3.1), we did not apply a redenning correction, but rather used the observed value, both here 
and in §4.2. 

■^Thc extended 12/im sample probes to a lower flux density limit than the original 12/im sample: 0.22 



Jy vs . 0.30 Jy, giving a total of 118 detected Sys over the 59 detected in the original sample (JRush et al. 



19931 ) 
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and [OIV] flux values are derived from Melendez et al. (2008) and Tommasin et al. (2008, 
2010), which presents high resolution Spitzer spectroscopy of the extended 12/ini sample. As 
only Winter et al. (2010) quote Balmer decrements, we only have comparison Syl F[oiii],corr 
values for the samples selected from the BAT catalog (i.e. Fi4_i95fcev and F[o/y] from Weaver 
et al. 2010). We utilize both the Kolmogorov-Smirnov test ("K-S" test) and Kuiper's test 
on the detected data points (excluding the three [OIV] and three radio upper limits in the 
12yum data) to determine to which extent the flux ratios are significantly inconsistent between 
the Syl and Sy2 galaxies: a lower Kuiper "D" value indicates that these two populations 
are drawn from the same parent population, suggesting that such fluxes are independent 
of viewing angle. The Kuiper test is similar to the more often used K-S test but with the 
following modification: the "D" statistic of the K-S test represents the maximum deviation 
between the cumulative distribution functions (CDFs) of the two samples, whereas the "D" 
statistic in Kuiper's test is the sum of the maximum and minimum deviations between the 
CDFs of the two samples, so that this statistic is as sensitive to the tails as to the median 
of the distribution. The results of the K-S test and Kuiper test agree in that they do 
not lead us to reject the null hypothesis that the two samples are drawn from the same 
parent population, with the exception of the F[oiv]/Fioiii],obs ratio, where the tests lead to 
conflicting results. We note that two-sample tests work better for larger data sets, so the 
probabilities quoted in Table [TOl should be interpreted as approximate. 

The comparisons of the Sy2 and Syl samples are shown in Figures IT] through [TOl In each, 
the dotted-dashed and (in the cases of more than one comparison sample) dashed line(s) on 
these plots indicate the mean values for the Syl diagnostic ratios and the correlations from 
linear regression we calculated from the literature values. 

A mild disagreement between the average Syl and Sy2 Fi^oiv]/F[oiii],obs ratio (up to 
a factor of ~2) is evident. Here the comparison Syl data come from the Diamond-Stanic 
et al. (2009) and the Melendez et al. (2008) samples. We obtain mixed results as to 
the significance of this difference based on the statistical test used: according to the K-S 
test, the F[o/i/]/F[o///],o6s ratio for the Syl and Sy2 populations are statistically significantly 
different (D=0.301, p=0.042), but not according to the Kuiper test (D=0.310, p=0.223). We 
find similar disparate results when we run these tests on the Fi^oiv]/F[oiii],obs ratio for the 
detected points between the Syls and Sy2s in the Melendez et al. (2008) and Diamond-Stanic 
et al (2009) samples, namely that the K-S test implies different parent populations (p=0.005 
and p=0.004, respectively) but not the Kuiper test (p=0.097 and p=0.126, respectively). 
Melendez et al. (2008) and Diamond-Stanic et al. (2009) (as well has Haas et al. 2005, who 
compared seven quasars with seven Fanaroff-Riley II (FRII) radio galaxies) have reported 
significant differences between the observed [OIII] and [OIV] flux between Syls (quasars for 
Haas et al. 2005) and Sy2s (FRIIs for Haas et al. 2005), with the type 2 sources having higher 
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'^ [oiv]/'^ [oiii],obs ratios. These authors have attributed the diminumtion of [OIII] in type 2 
AGN to extinction in the NLR. Baum et aL (2010) suggests that such [OIII] obscuration 
results from the AGN torus: using the 12 /im sample, they find a correlation between the 
F[o/y]/ ^[oni],obs ratio and the Sil 10/xm feature, which probes torus obscurationjj In type 
1 Syls, this silicate feature is in emission, whereas Sy2s exhibit Sil absorption, making the 
Sil strength a probe of system orientation. The ratio of ^[oiv] to ^[oiii],ohs increases with Sil 
absorption (parameterized by negative values of the Sil strength) which could suggest that 
the torus is extincting part of the [OIII] emission. Our results may confirm these previous 
studies as we find that Sy2s tend to have lower observed [OIII] emission as compared to Syls 
and this may be due to NLR extinction. We note, however, that such extinction affects the 
[OIII] line only up to a factor of 2 on average between our Sy2 and comparison Syl samples, 
albeit with a wide dispersion, and this difference between the two populations may not be 
significant. 

The average log (P MiR/^[oni],obs) ratio is consistent between Syls (2.56 ± 0.50 dex) and 
Sy2s (2.62 ± 0.62 dex), which could seemingly contradict the results cited above where NLR 
extinction causes attenuation of the [OIII] flux in Sy2s but not in Syls. The clumpy torus 
model of Nenkova et al. (2008) and smooth torus model of Pier & Krolik (1992) predicts 
a slight anistropy in emission at 12/im depending on viewing angle: as the viewing angle 
increases from 0° (Syl) to 90° (Sy2), the torus flux decreases by a factor of ~2. The effects of 
depressed MIR emission in Sy2s and enhanced MIR emission in Syls, assuming [OIII] is more 
extincted in the former than the latter, would therefore result in 'PMiR/'P[oiii],ohs ratios that 
are more consistent than 'P[oiv]/'P[oiii],ohsi which is indeed what we observe. However, the 
average differences between PMiR/P[oiii],obs and P[oiv]/P[oiii],obs are within the scatter of 
these ratio values, and we are unable to rule this out as the main driver for the disagreement, 
rather than invoking anisotropies in torus emission. 

Interestingly, the relationship between Ijmir and Ljo/y] is nearly identical for Syl and 
Sy2 galaxies (Figure [3]). Though the MIR flux is not corrected for starburst contamination 
(see Appendix), and the Sy2s in the 12/xm sample are thought to harbor more star formation 
activity than Syls (e.g. Buchanan et al. 2006), we see no evidence that star formation 
activity is contributing significantly to the MIR emission. As the 'Pmir/'P[oiv] diagnostic 
ratio shows the smallest dispersion for the combined sample, a similar relationship to Syl 
galaxies and no evidence for luminosity bias (see Sectin §4.3), the MIR and [OIV] flux may 
be the most robust proxies for the intrinsic AGN luminosity in Type 2 AGN. However, the 
KS test and Kuiper's test indicates a lower probability that the Syl and Sy2 samples agree 



^They define the Sil strength by the natural logarithm of the observed flux of the feature divided by the 
intcropolatcd continuum flux at 10/im. 
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than the 'PMiR/^[oiii],obs ratio, though this could be driven by the lower Ng// value for 
Fm/_r/F[o///],o6s rather than a robust statistical agreement. 

The different slopes between Syls and Sy2s in the luminosity plots of the radio data 
against other intrinsic AGN flux proxies (Figures HI |5] and [6]) suggest disagreements between 
these samples. However, the 'P[oiv]/'P%aghz and ^mir/^%aghz flux ratios are consistent 
between Syl and Sy2 galaxies, indicating that the disparate slopes are perhaps influenced 
by scatter due to the wide range of radio loudness in AGN. Results of the KS test and 
Kuiper's test (Table [TU1) also indicate that the differences in the radio flux ratios between 
Syls and Sy2s are not statistically significant. Diamond-Stanic et al. (2009) compared the 
6 cm radio data between Syls and Sy2s and found that for the Sy2s with a measured X- 
ray column density, these two samples show no statistically significantly differences, though 
they find a higher probability that they are drawn from the same distribution (~68 - 78%) 
than we do (~55%). Melendez et al. (2010) also found that the 8.4 GHz and [OIV] fluxes 
between Syl and Sy2 galaxies are not significantly different, though sources dominated by 
star formation (i.e. less than 50% of the [Nell] line attributable to AGN ionization) had 
statistically different ^[oiv]/^%aghz values than AGN dominated sources, indicating that 
radio emission may not accurately trace intrinsic AGN power. This latter result may agree 
qualitatively with our Figure |5l where the "PAH-strong" sources lie at or below the best-fit 
line between L[o/v'] and \jsaghz- 

The hard X-ray proxy performs much more poorly (Figures 7 - [TUl) . based on both the 
wider dispersion in the diagnostic flux ratios and the larger disagreement between the Syl 
and Sy2 flux ratios. The mean hard X-ray emission (normalized by other isotropic indicators) 
in Sy2s tends to be about an order of magnitude weaker than in Syls, though this is driven 
largely by the 12/im sample as only one source was detected by BAT in the [OIII] sample. 
This disagreement agrees with the results of Rigby et al. (2008) and Weaver et al. (2010), 
where the X-ray flux was normalized by the [OIV] emission. Indeed, using survival analysis, 
we flnd that Fi4_i95fcey/F[o/y] disagrees signiflcantly between BAT-selected Syls and both 
the [OIII] and 12/im sub-samples. Such a large disagreement is not found between the Syls 
and Sy2s in the Winter et al. (2010) sample (see Table [9]), which is driven by the lower [OIII] 
flux observed in their Sy2s as compared to their Syls. The hard X-ray to [OIII] flux ratios, 
both observed and reddening corrected, do not differ signiflcantly between the BAT-selected 
Syls and the [OIII]-selected Sy2s, but do for the 12/im sample[] According to the Logrank 



^Fi4_i95fcey/F[o7//]_o&s has mixed results for the 12/im sample: according to two of the tests (Logrank 
and Peto & Prentice Generalized Wilcoxon Test), the difference is significant though with other tests, the 
null hypothesis that the two samples are drawn from the same parent distribution can only be discarded at 
the ~12% confidence level. 
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and Peto & Prentice Generalized Wilcoxon tests, the hard X-ray flux normalized by the MIR 
flux differs signficantly for both Sy2 subsamples and the BAT-selected Syls. Consistent 
with the results from §4.1, hard X-ray selected AGN do not represent the population of 
those selected in the MIR, and there may be some evidence that they do not fully sample 
the optically selected sources. Compton scattering may be responsible for weakening the 
observed hard X-ray emission in Sy2s, as suggested by Weaver et al. (2010), which indicates 
that the 14-195 keV emission is not truly isotropic. 



4.3. Luminosity Dependence 

As we have seen above, there is significant scatter in the flux ratios of the different 
proxies for AGN intrinsic luminosity. Here we examine the possibility that some of this 
scatter is caused by systematic differences that correlate with the accretion rate of the black 
hole (in units of the Eddington limit). 

To make this test, for any pair of luminosity proxies we parameterized Lagn /LEdd by 
the square root of the product of the luminosities of the two proxies divided by the mass of 
the central black hole {M-bh-, listed in Tables [1] and |2]) . Linear regression fits were performed, 
with the correlation coefficients and probability of uncorrelation listed in Table [TTJ 

We find three statistically significant anti-correlations (Figure [13]): ^[oiv]/^[oiii],obs, 
FMrn/Fioiiiiobs and Fmir/Fsaghz- The anti-correlations for the ratios involving Fioiii],obs 
are largely driven by those galaxies with a high Balmer decrement. When we exclude the 6 
sources with Ha/H/3 > 9, which may be those systems with the most NLR etinction, these 
anti-correlations are no longer statistically significant. This may indicate that the bolometric 
correction to the observed [OIII] luminosity might have a weak dependence on the Eddington 
ratio. However, the observed [OIII] luminosity, which partly parameterizes the Eddington 
ratio, does not as accurately trace intrinsic AGN flux for these more dust-obscured sources. 
If the Eddington ratio is defined as just Lqiv I^bh and Lmir/Mbh in these relationships, 
the anti-correlations are no longer significant. Hence, the weak trends in Figure [T^ a) and b) 
is likely driven more by NLR extinction bias on the [OIII] flux rather than the accretion rate 
of the black hole. This latter mechanisms could be affecting the Fmir/Fsagh ratio, albeit 
with wide scatter. 
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5. Starburst Activity in the Host Galaxy 

5.1. Comparison of Different Diagnostics 

Given the strong connection between Type 2 AGN and star-formation (e.g. Kauffmann 
et al. 2003) we expect that the signature of both processes will be present in optical and MIR 
spectra of AGN. By analogy to the previous section (where we compared different proxies 
for the intrinsic AGN luminosity) we now undertake a comparison of different diagnostics of 
the relative energetic significance of a starburst vs. the AGN. 

One such diagnostic involves the use of the MIR polycyclic aromatic hydrocarbon (PAH) 
features. These have been shown to be correlated with star formation activity (e.g. Smith 
et al. 2007) and possibly anti-correlated with the presence of an AGN (O'Dowd et al. 2009; 
Voit 1992). More specifically, we used the equivalent width (EW) of the PAH features to 
assess the relative amount of starburst activity in the host galaxy (e.g. Genzel et al. 1998). 
Another empirical diagnostic of the relative contribution of the starburst in the MIR can 
be parametrized by the MIR spectral index: a2o-30AtJj- Larger values of 0^20-30/^™ indicate 
the presence of cold dust from starburst activity (Deo et al. 2009 and references therein). 
The ratio of the [OIV] to [Nell] 12.81/im MIR emission-lines probes the hardness of ionizing 
spectrum and hence the relative importance of the AGN and starburst. A larger ratio (~1) 
implies t he dominance of A GN activity whereas a lower ratio (~0.02) is pure starburst 



activity (jGenzel et al. Ill998l ). The analogous diagnostic from optical spectra is the distance 
a galaxy spectrum lies from the locus of star forming galaxies in the Baldwin, Phillips 
& Televich BPT (1981, BPT) diagram (D = ^{[NII]/Ha + 0.45)^ + {[0III]/Hl3 + 0.5)^, 
Kauffmann et al. 2003). A larger "D" parameter indicates pure AGN activity while a smaller 
value implies a mixture of starburst and AGN processes in the host galaxy. 

Figures [T3] and [T5] illustrate the relationship between these AGN and star formation 
activity indicators for the Sy2s in our combined sample. The color coding is the same as 
in Figures [1] - [TOl According to linear regression analysis, a2o-3o^m and the PAH EWs are 
correlated and [OIV]/[NeII] and the PAH EWs are anti-correlated at greater than the 3.5cr 
level: PAH EW 11.3 /im vs. a2o^3o^^m has p=0.609, P„„corr=l-47xl0-5; PAH EW 17 /im 
vs. a2o-3o^.m has p=0.600, Puncorr=5.26xl0~^; PAH EW 11.3mMm vs. [OIV]/[NeII] has p=- 
0.677, Funcorr=2.26xl0~^; PAHEW 17/imvs. [OIV]/[NeII] hasp=-0.515, P„„corr=4.89xl0-^ 
We also note that the Sy2s with strong PAH emission mostly lie at systematically higher 
PAH EW values than the relation found between a2o-30fim and the PAH EWs. The majority 
of the Sy2s have D values between ~1. 2-2.0, with five of the strong PAH sources at system- 



8 



aAi-Ao = log(/Ai//A2)/log(Ai/A2) 
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atically lower D values, ~0. 5-1.1 (see Figure [Tarl ). The Sy2s with lower D values have higher 
PAH EW values and lower 'P[oiv]/'^[Neii] ratios, though they exhibit a weaker trend with 
the IR spectral index. 

These results indicate that the various MIR and optical indicators of starburst activity 
agree both qualitatively and quantitatively. 



5.2. The Spatial Scale of the MIR Emission and the Role of the Host Galaxy 

The results above refer to measurements of the central region enclosed by the IRS 
aperture. However, the 12/im sample was drawn from the IRAS Point Source Catalog, where 
the aperture size (0.75 x 4.5' at 12/im) is much larger and will encompass contributions from 
the host galaxy. To quantify the extendedness of the MIR emission in the 12/im sample, we 
calculate the ratio of the IRAS flux (at 12/tm, Spinoglio & Malkan 1989) to the IRS flux. A 
ratio of ~1 indicates the MIR emission is dominated by the galactic center, whereas higher 
ratios imply a greater amount of extended emission. In Figure [T71 we plot this ratio against 
the PAH EW at 11.3 /im and 'P[oiv]/'P[Neii]- As expected, the relative amount of extended 
MIR emission decreases as the relative energetic importance of the AGN increases. 

Comparison of the ratio by which the SL module was rescaled to match the LL module 
(see § 4.2) reveals the presence of extended emission on smaller spatial scales. Inspection of 
this extendedness factor does not show any signiflcant differences between the "PAH-weak" 
and "PAH-strong" sources (with the exception of NGC 7582). 



5.3. Are [OIII] and [OIV] Biased by Star Formation? 

In this section, we investigate whether the relative fraction of the AGN bolometric 
luminosity that emerges in [OIII] and [OIV] line emission is preferentially higher in galaxies 
with more star formation. This might be expected if the gas clouds in the NLR that are 
photoionized by the AGN and produce these lines are directly related to the gas clouds 
responsible for star-formation. If true, this would imply that AGN selected using [OIII] or 
[OIV] would be biased towards galaxies with higher star formation rates. 

To test this, we have plotted the ratio of both ^[oiii],ohs and 'P[oiv] to Fm/r versus 



^Thc two "PAH-weak" sources with low D values are Arp 220 and F08572+3915; which have [OIV] upper 
limits and no measureable PAH EW at 11.3^m. 
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the star formation indicators analyzed above (PAH EWs, IR spectral index and the "D" 
parameter). We find no strong trends between the star formation indicators and [OIV] and 
[OIII] emission, as illustrated in Figures [TS] and ^U\ 

We conclude that there is no convincing evidence that host galaxies with a large star 
formation rate have preferentially higher relative luminosities of [OIII] and [OIV] at the 
luminosities represented in this sample, where the bolometric luminosity {l-n,oi) ranges from 
Uoi ?« 10^ - 8 X 10^^ Lsun, which is ~3xl0~^ to 0.5 of the Eddington luminosity (Leddfl 
Thus, these proxies of intrinsic AGN power are not biased by star formation activity at these 
Eddington ratios. 



6. Conclusions 

We have taken an empirical approach in analyzing the agreement among the various 
indicators of isotropic AGN luminosity for two complete and homogeneously selected sam- 
ples of Sy2s, one selected based on observed [OIII] flux and the other on MIR flux. The 
diagnostic ratios with the smallest spread are likely those where such biases from sample 
selection, starburst contamination, statistical errors, and scatter due to the various physical 
mechanisms that produce these emission features, are minimized. Such indicators, as well 
as those that agree the most with Syl relations, may be the most robust tracers of AGN 
activity. Our results on these indicators are summarized below. 

- Sample Selection The optically selected sample, picked on the basis of high [OIII] flux, 
shows tighter correlations among these diagnostics than the MIR selected sample. We inves- 
tigated whether the inclusion of active star forming galaxies in the 12/im sample introduced 
the spread in these ratios by dividing the sample into galaxies that have large amounts of 
starburst activity ( "PAH-strong" ) and those that do not ( "PAH- weak" ) . The distribution 
of the diagnostic ratios for the two sub-samples are similar. Isolating the 12/im sources with 
a limited D range (1.2<D<1.7) also results in large scatter that is systematically higher 
than that observed in the [OIII] sample for all but three of the flux ratios. A similarly wide 
spread between Fioiv]/F[oiii],obs is present in other samples (i.e. Melendez et al. 2008 and 
Diamond-Stanic et al. 2009), indicating that sample selection based on [OIII] is primarily 
responsible for the tighter correlations we observe. This may be due to less extinction in the 
NLR region which would be expected in sources that have high observed [OIII] flux. 



i^We estimated Lfcp; by a ssuming the observed mid-infrared flux constitutes 20% of the bolometric flux 
dSpinoglio fc Malkan iflisol) . 
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- Extinction Correction Applying an extinction correction to the [OIII] flux tightens the 
correlations with the other luminosity indicators for the [OIII]-selected sample, yet broadens 
the dispersion for the 12/im sample. It is not clear whether this difference is primarily due to 
the different sources of the emission-line data (homogenous SDSS data for the [OIII] sample 
and heterogeneous data for the 12/im sample), or whether it simply points out the limitations 
of extinction corrections in very dusty AGN. Comparison of the optical vs. MIR properties 
of the most dusty AGN in the SDSS suggest that the former effect is important (Wild et al. 
2010). 

- Agreement Among Sy2s The observed [OIII], [OIV] and MIR luminosities agree the 
best in the combined Sy2 sample. The widest spread among the various proxies is seen in the 
radio emission. The X-ray data are dominated by upper limits, but also show a significantly 
larger dispersion than the optical and IR isotropic flux indicators. 

- Comparison with Syls The mean ratio of the observed [OIII] flux to the [OIV] flux is 
lower in Sy2 than in Syls by a factor of 1.5-2, while the mean ratio of the observed [OIII] 
flux to MIR is consistent between Sy2s and Syls. The former result, which represents a 
statistically significant difference between Syls and Sy2s according to the KS test but not 
Kuiper's test, agrees with previous findings (e.g. Haas et al. 2005, Melendez et al. 2008, 
Diamond-Stanic et al. 2009) and has been interpreted as extinction affecting [OIII] in the 
NLR, or even torus obscuration attenuating the [OIII] emission (Baum et al. 2010). However, 
the latter result can not be simply explained by a larger amount of dust extinction of [OIII] in 
the Sy2s, but it could be due to a slight anisotropy in the MIR emission as predicted by Pier & 
Krolik (1992) and Nenkova et al. (2008) where Syls could have up to a factor of two higher 
MIR flux as compared to Sy2s. The wide scatter in these ratios can also be responsible 
for the discrepancy between the F[o/y]/F[o///],ofes and FMiB/F[oiii],obs values, which may 
be the main driver for the mild disagreement rather than torus emission anisotropy. The 
F[oiv]/F8.4:GHz and Fmir/F8 4^ghz mean values are consistent between Syl and Sy2 galaxies 
(in agreement with Diamond-Stanic et al. 2009 and Melendez et al. 2010 for the [OIV] and 
radio comparison), though the slopes of the luminosity plots show disagreements and there 
is wide scatter which is expected due to the wide range of radio loudness observed in AGNs. 

- Hard X-ray Selected Samples We find that the hard X-ray flux (relative to the [OIII], 
[OIV], and MIR fluxes) is suppressed by about an order-of-magnitude in MIR selected Sy2s 
compared to both Syls (consistent with Rigby et al. 2008) and to hard X-ray selected Type 
2 AGN (Winter et al. 2010 and Weaver et al. 2010). The comparison with the [OIII] sample 
is mixed, with statistically significant differences between the Sy2s and Syls when the X- 
ray flux is normalized by the [OIV] and MIR flux, but not when normalized by the [OIII] 
flux. However, hard X-ray selected Sy2s differ signiflcantly from [OIII]-selected Sy2s (though 
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with only one [OIII] Sy2 detected by BAT, this analysis may be less robust than the 12/im 
comparison). These results indicate that hard X-ray emission (E > 10 keV) is anistropic and 
hard X-ray selected samples are biased against the more heavily obscured type 2 AGN that 
are present in MIR and possibly [OIII] selected samples. As Weaver et al. (2010) note for 
sources detected in hard X-rays, Compton scattering could be responsible for the hard X-ray 
attenuation observed in Type 2 AGN as compared to Type 1. In more obscured sources, 
Compton scattering may be pushing them below the flux sensitivity of BAT. 

Fmir and F[o/y] agree the best, both in comparison with the other indicators in the 
combined Sy2 sample (having the least scatter) and in having a nearly identical flux ratio 
in Sy2s and Syls as well as not suffering from a luminosity bias. Among the indicators we 
have considered, they are the best proxies for truly isotropic AGN emission. 

We tested the level of agreement of various optical and infrared indicators of star for- 
mation activity compared with proxies of AGN activity. Similar to previous works, we find 
statistically significant relations among the various indicators of the relative energetic signif- 
icance of star formation and an AGN. These include the MIR spectral slope (parametrized 
by «2o-30Mm)! the PAH EWs at 11.3yU.m and 17yU.m, the ratio of [OIV]/[NeII] fluxes, and the 
location of the galaxy in the commonly used diagnostic diagram based on optical emission- 
lines. We note that the last two diagnostics are a measure of the relative contribution of 
AGN vs. starburst activity to the incident ionizing radiation. As the incident radiation field 
becomes more dominated by the AGN activity, the PAHs become weaker relative to the 
MIR continuum. In part this is expected because of the "dilution" of the MIR continuum 
by AGN-heated dust, but the AGN may also be directly affecting the PAH emission (e.g. 
O'Dowd et al. 2009). 

We also found that the Sy2s that were clearly starburst/AGN composites based on the 
above indicators were systematically those cases in which large-scale MIR emission from the 
host galaxy greatly exceeded that from the AGN. We quantified this by comparing the ratio 
of the 12/im flux from the large IRAS aperture with the 15.5/im flux from the small IRS 
aperture. A smaller aperture is therefore necessary to isolate AGN emission in systems with 
high rates of star formation on large scales and/or low AGN luminosities. 

The ratios of the [OIII] /MIR and [OIV]/MIR fluxes show little if any evidence for a 
correlation with any of the above measures of the relative amount of star formation. This lack 
of a relationship suggests that the [OIII] and [OIV] lines are not biased to be a preferentially 
higher fraction of the AGN bolometric luminosity in host galaxies with more star formation 
activity (more dense gas). 

This work is based in part on observations made with the Spitzer Space Telescope, which 
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is operated by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. The authors thank the anonymous referee whose insightful comments 
improved the quality of this manuscript. 
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Fig. 1. — a) Distribution of log {P[oiv]/'P[oiii],corr)- Black solid line shows the combined 
sample, filled cyan histogram depicts the [OIII]-sample, and red (green) dashed histogram are 
the weak PAH (strong PAH) 12/im sources, b) Log {hpiv]) vs log {h[oiii],corr)] black dashed 
line represents best-fit from linear regression with p=0.837 giving P„„corr=l-28xl0~^^ (slope 
= 0.83 ± 0.08 dex with intercept at 6.32 dex). Cyan asterisks represent the [OIII] sample 
and red diamonds (green triangles) depict the weak PAH (strong PAH) 12/xm Sy2s. The blue 
dotted-dashed line in a) and b) refiect the values for Syls, with 'P[oiii],corr from Winter et al. 
(2010) and F[o/y] from Weaver et al. (2010). c) Distribution of log {P [oiv] I'P [oiii],ohs) ■ d) 
Log(L[o/\/]) vs log(L[o///],ofes); black dashed line represents best-fit from linear regression with 
p=0.884 giving P„„,corr=8.52x 10"^^ (slope = 0.86 ± 0.07 dex with intercept at 5.92 dex). The 
purple dashed line and blue dot-dashed line represent the Syl values from Diamond-Stanic 
et al. (2009) and Melendez et al. (2008) samples, respectively. 
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Fig. 2.— a) Distribution of log (FMiR/F[oiii],corr)- b) Log (Lmir) vs log moiii],corr)] 
p=0.711 giving P„„corr=l-04xl0"® (slope = 0.71 ± 0.10 dex with intercept at 13.8 dex). 
c) Distribution of log (FMiR/F[oiii],obs)- d) Log (Lmir) vs log {Lioiii],obs); p=0.731 giving 
P«ncorr=2.46xl0~^ (slope = 0.68 ± 0.09 dex with intercept at 15.7 dex). The blue dashed- 
dotted line represents the Syl values, with Fmir from Deo et al. (2009) and Fi^oiii],obs from 
Melendez et al. (2008). Color and linestyle coding same as Figured! 
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and Ffo/y] from Tommasin et al. (2010) and Melendez et al. (2008). The relationship 
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same as Figure [TJ 
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Fig. 4.— a) Distribution of log (F[oiii],corr/FsAGHz)- b) Log (L[o///],corr) vs log (L8,4g//2); 
p=0.667 giving P„„corr=2.00xl0~^ (slope = 0.57 ± 0.13 dex with intercept at 20.1 dex). c) 
Distribution of log {F[oiii],obs/FsAGHz)- d) Log moiii],obs) vs log (Ls^^ghz)] p=0.57 giving 
Puncorr=0.003 (slope = 0.49 ± 0.14 dex with intercept at 22.1 dex). The blue dashed-dotted 
line represents the Syl values, with the radio data from Thean et al. (2000) and F[oiii],obs 
from Diamond-Stanic et al. (2009) and Melendez et al. (2008). Color and linestyle coding 
same as Figure [H 
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Table 1: [OIII] Sample 



Galaxy 



RA Dec 

(J2000) (J2000) 



Distance 

Mpc 



Mo 



NGC 0291 


00 53 29.9 


08 46 04 


0.019 


82.6 


7.10 


Mrk 0609 


03 25 25.4 


06 08 37 


0.034 


151.3 


7.80 


IC 0486 


08 00 21.0 


26 36 49 


0.027 


116.8 


7.73 


2MASX J08035923+2345201 


08 03 59.2 


23 45 20 


0.029 


128.8 


7.25 


2MASX J08244333+2959238 


08 24 43.3 


29 59 24 


0.025 


111.0 


6.99 


CGCG 064-017 


09 59 14.8 


12 59 16 


0.034 


150.9 


6.77 


2MASX J10181928+3722419 


10 18 19.3 


37 22 42 


0.049 


219.8 


7.14 


2MASX J11110693+0228477 


11 11 06.9 


02 28 48 


0.035 


154.9 


7.33 


CGCG 242-028 


11 23 01.3 


47 03 09 


0.025 


110.1 


7.29 


SBS 1133+572 


11 35 49.1 


56 57 08 


0.051 


229.0 


7.99 


Mrk 1457 


11 47 21.6 


52 26 58 


0.049 


216.6 


7.85 


2MASX J11570483+5249036 


11 57 04.8 


52 49 04 


0.036 


156.7 


7.17 


2MASX J12183945+4706275 


12 18 39.4 


47 06 28 


0.094 


431.4 


8.17 


2MASX J12384342+0927362 


12 38 43.4 


09 27 37 


0.083 


377.9 


8.40 


CGCG 218-007 


13 23 48.5 


43 18 04 


0.027 


119.4 


7.65 


2MASX J13463217+6423247 


13 46 32.1 


64 23 24 


0.024 


105.6 


6.81 


NGC 5695 


14 37 22.1 


36 34 04 


0.014 


61.1 


7.44 


2MASX J15522564+2753435 


15 52 25.7 


27 53 43 


0.075 


337.6 


8.05 


SBS 1609+490 


16 10 51.8 


48 54 39 


0.045 


200.0 


6.87 


2MASX J16164729+3716209 


16 16 47.3 


37 16 21 


0.152 


725.0 


8.34 



Distances based on redshift, assuming a cosmology of 
Hq = 70 km s-\ Qm = 0.27 and Q^ac = 0.73. 
^Masses based on velocity dispersion: 
Mbh = 10^-^'^{a/200 km s-^)^-°2 M© fTremaine et al. 



2002). 
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Table 2: 12;um Sample 



Galaxy 


RA 


Dec 


z 


Distance 


Mbh' 


Ref2 




(J2000) 


(J2000) 




Mpc 


Mq 




IRAS 00198-7926 


00 21 53.6 


-79 10 07 


0.073 


325.8 


8.44 


B07 


NGC 0424 


01 11 27.6 


-38 05 00 


0.012 


51.2 


7.78 


B07 


NGC 1068 


02 42 40.7 


-00 00 48 


0.004 


16.9 


7.35 


B07 


NGC 1143/4 


02 55 09.7 


-00 10 40 


0.028 


120.8 


8.29 


B07 


NGC 1320 


03 24 48.7 


-03 02 32 


0.009 


38.3 


7.18 


B07 


NGC 1386 


03 36 46.2 


-35 59 57 


0.003 


12.7 


7.24 


B07 


F04385-0828 


04 40 54.9 


-08 22 22 


0.015 


64.1 


8.77 


W07 


NGC 1667 


04 48 37.1 


-06 19 12 


0.015 


64.1 


7.88 


B07 


F05189-2524 


05 21 01.5 


-25 21 45 


0.043 


187.7 


7.86 


W07 


F08572+3915 


09 00 25.4 


39 03 54 


0.058 


256.0 


7.95 


V09 


NGC 3982 


11 56 28.1 


55 07 31 


0.004 


16.9 


6.09 


B07 


NGC 4388 


12 25 46.7 


12 39 44 


0.008 


34.0 


7.22 


B07 


NGC 4501 


12 31 59.2 


14 25 14 


0.008 


34.0 


7.86 


B07 


TOLOLO 1238-364 


12 40 52.8 


-36 45 21 


0.011 


46.9 


6.84 


W07 


NGC 4968 


13 07 06.0 


-23 40 37 


0.010 


42.6 


7.27 


P06 


MCG-3-34-64 


13 22 24.4 


-16 43 43 


0.017 


72.7 


7.70 


P06 


NGC 5135 


13 25 44.0 


-29 50 01 


0.014 


59.8 


7.35 


B07 


NGC 5194 


13 29 52.7 


47 11 43 


0.002 


8.5 


6.95 


B07 


NGC 5347 


13 53 17.8 


33 29 27 


0.008 


34.0 


6.79 


B07 


Mrk 463 


13 56 02.9 


18 22 19 


0.050 


219.4 


7.88 


B07 


NGC 5506 


14 13 14.9 


-03 12 27 


0.006 


25.5 


6.65 


B07 


NGC 5929 


15 26 06.1 


41 40 14 


0.008 


34.0 


7.25 


B07 


NGC 5953 


15 34 32.4 


15 11 38 


0.007 


29.7 


9.80 


P06 


Arp 220 


15 34 57.1 


23 30 11 


0.018 


77.1 


7.78 


V09 


NGC 6890 


20 18 18.1 


-44 48 25 


0.008 


34.0 


6.48 


B07 


IC 5063 


20 52 02.3 


-57 04 08 


0.011 


46.9 


7.74 


B07 


NCG 7130 


21 48 19.5 


-34 57 05 


0.016 


68.4 


7.53 


B07 


NGC 7172 


22 02 01.9 


-31 52 11 


0.009 


38.3 


7.67 


B07 


NGC 7582 


23 18 23.5 


-42 22 14 


0.005 


21.2 


7.25 


B07 


NGC 7590 


23 18 54.8 


-42 14 21 


0.005 


21.2 


6.79 


B07 


NGC 7674 


23 27 56.7 


08 46 45 


0.029 


125.3 


7.56 


B07 



Distances based on redshift, assuming a cosmology of 

Hq = 70 km s-\ Q.M = 0.27 and 9.^ac = 0.73. 

Coordinates and redshift values from NED. 

^Masses based on velocity dispersion: Mbh = 10^-^^(o-/200 km s"^)^-^^ M© fTremaine et al. 2002) . 

except for F08572 -1-3915 which is based on photometry of the host galaxy (see V09 for details). 

The FWHM of the [OIII] line was used as a proxy for a for F04385-0828, 

F05189-2524 and TOLOLO 1238-364 (Wang & Zhang 2007, Greene & Ho 2005). 

^References: B07 - Bian & Gu 2007; P06 - Peng et al. 2006; 

V09 - Veilleux et al. 2009; W07 - Wang k Zhang 2007. 
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Table 3: Optical Emission Line Fluxes(10 erg s cm ) and Ratios for [OIII] Sample 



Galaxy 


[OIII] Ob, 


[Olll]carr 


Ha/H/3 


[Nil] /Ha 


[0III]/H/3 


NCG 0291 


6.47 


30.4 


5.06 


0.70 


7.94 


Mrk 0609 


7.54 


31.0 


4.85 


0.93 


5.71 


IC 0486 


4.59 


17.1 


4.70 


1.17 


9.34 


2MASX J08035923+2345201 


5.22 


6.45 


3.32 


1.02 


10.1 


2MASX J08244333+2959238 


8.76 


19.31 


3.98 


0.37 


9.30 


CGCG 064-017 


8.46 


12.33 


3.49 


0.24 


6.63 


2MASX J10181928+3722419 


4.05 


6.01 


3.51 


0.84 


11.2 


2MASX J11110693+0228477 


4.04 


5.45 


3.41 


0.91 


9.31 


CGCG 242-028 


4.50 


7.92 


3.71 


0.75 


7.66 


SBS 1133+572 


4.71 


27.3 


5.41 


0.90 


6.28 


Mrk 1457 


4.27 


16.2 


4.72 


0.78 


4.63 


2MASX Jl 1570483+5249036 


4.95 


7.17 


3.49 


0.54 


11.1 


2MASX J12183945+4706275 


4.94 


10.4 


3.93 


0.21 


10.6 


2MASX J12384342+0927362 


5.48 


10.2 


3.77 


0.50 


11.3 


CGCG 218-007 


4.12 


21.3 


5.22 


0.96 


6.24 


2MASX J13463217+6423247 


4.15 


9.60 


4.04 


1.24 


10.1 


NGC 5695 


5.93 


9.30 


3.57 


1.36 


9.56 


2MASX J15522564+2753435 


5.87 


7.69 


3.38 


0.57 


12.4 


SBS 1609+490 


4.68 


8.22 


3.71 


0.86 


10.6 


2MASX J16164729+3716209 


5.03 


7.21 


3.47 


0.35 


10.3 



Table 4. Optical Emissiori 


L Line Fluxes (10"^^ 


erg s ^ cm ^) and Ratios for 


12/im Sampl 


e 




Object 


[OIII]o,, 


yJlLl corr 


Ha/H/3 


Refi 


Nil] /Ha 


Ref 


0III]/H/3 


Ref 




IRAS 00198-7926 


3.88 


38.5 


6.41 


D92 


0.63 


V93 


11.5 


V93 




NGC 0424 


70.0 


121 


3.69 


V97 


0.30 


V97 


4.55 


V97 




NGC 1068 


1367 


21280 


7.39 


S95 


1.80 


S95 


13.0 


S95 




NGC 1144 


4.30 


28.7 


5.65 


D88 


1.73 


M06 


7.10 


D88 




NGC 1320 


14.0 


58.0 


4.86 


T03 


0.60 


D92 


8.98 


D92 




NGC 1386 


79.2 


1116 


7.16 


S89 


1.28 


S89 


16.1 


S89 




F04385-0828 


0.35 


10.9 


9.20 


T03 


>0.65 


D92 


>2.15 


D92 




NGC 1667 


6.38 


238 


9.74 


S95 


0.69 


S95 


3.48 


S95 




F05189-2524 


15.3 


348 


8.33 


V99 


0.46 


V99 


44.0 


V99 




F08572+3915 


0.12 


1.41 


6.71 


SDSS 


0.40 


SDSS 


0.97 


SDSS 


1 


NGC 3982 


20.0 


64.9 


4.50 


T03 






17.4 


T03 


1 


NGC 4388 


80.0 


557 


5.73 


D92 


2.47 


D92 


11.8 


D92 




NGC 4501 


3.70 


5.48 


3.51 


T03 


1.86 


S82 


5.82 


S82 




TOLOLO 1238-364 


45.7 


289 


5.56 


S95 


0.69 


S95 


5.31 


S95 




NGC 4968 


17.7 


1255 


11.9 


D92 


1.01 


D92 


21.5 


D92 




M-3-34-64 


151 


3190 


8.14 


D92 


0.99 


D92 


19.4 


D92 




NGC 5135 


37.5 


738 


7.96 


S95 


0.87 


S95 


4.30 


S95 




NGC 5194 


5.14 


45.0 


6.16 


M06 


2.40 


M06 


6.43 


M06 




NGC 5347 


4.43 


16.2 


4.67 


SDSS 


0.66 


SDSS 


7.79 


SDSS 




Mrk 463 


56.3 


150.3 


4.23 


D92 


0.40 


D92 


9.01 


D92 




NGC 5506 


127 


959 


5.88 


WIO 


0.70 


WIO 


7.29 


WIO 




NGC 5929 


6.92 


21.6 


4.44 


SDSS 


0.59 


SDSS 


3.42 


SDSS 




NGC 5953 


4.73 


80.1 


7.59 


V95 


0.78 


V95 


2.00 


V95 




Arp 220 


0.18 


33.7 


16.2 


SDSS 


1.41 


SDSS 


2.33 


SDSS 




NGC 6890 


23.3 


86.1 


4.69 


S90 


0.76 


S90 


10.2 


S90 











Table 4— 


Continued 








Object 


[0III]o5s 


[OIII],o,, 


Ha/H/3 


Refi 


Nil] /Ha 


Ref 


0III]/H/3 


Ref 


15063 


175 


224 


3.35 


S90 


0.54 


S90 


9.71 


S90 


NGC 7130 


34.1 


668 


7.95 


S95 


0.75 


S95 


4.84 


S95 


NGC 7172 


3.70 


214 


11.2 


V97 


0.30 


V97 


10.0 


V97 


NGC 7582 


51.3 


385 


5.87 


S95 


0.78 


S95 


3.21 


S95 


NGC 7590 


2.91 


17.3 


5.45 


S95 


0.65 


S95 


1.07 


S95 


NGC 7674 


71.8 


366 


5.19 


D92 


0.64 


D92 


11.3 


D92 






^Optical References: 
D88: Dahari & De Robertis, 1988; D92: de Grijp et al. 1992; M06: Moustakas & Kennicutt 
2006; S89: Storchi-Bergmann et al. 1989; S90: Storchi-Bergman et al. 1990; S95: Storchi- 
Bergmann et al. 1995; T03: Tran 2003; V97: Vaceli et al 1997, V95: Veilleux et al 1995; 
V99: Veilleux et al. 1999; WIO: Winter et al. 2010 
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Table 5. Spitzer Observation Summary^ 



Galaxy 


Program ID 


Mode 


Extraction Area^ 




5) 


kpc 


IRAS 00198-7926 


30291^ 


Staring 






NGC 0424 


30291^ 


Staring 








03269 


Mapping 


15.3x20.4 


3.5x4.6 


NGC 1068 


00014 


Staring 






NGC 1144 


00021^ 


Staring 








00021 


Mapping 


20.4x20.4 


11.3x11.3 


NGC 1320 


30291 


Staring 








03269 


Mapping 


15.3x20.4 


2.6x3.5 


NGC 1386 


00086 
30572^ 


Staring 
Staring 






F04385-0828 


30291^ 


Staring 








03269 


Mapping 


15.3x20.4 


4.7x6.2 


NGC 1667 


30291^ 


Staring 








03269 


Mapping 


15.3x20.4 


4.7x6.3 


F05189-2524 


00105 


Staring 






F08572+3915 


00105 


Staring 






NGC 3982 


30291^ 


Staring 








03269 


Mapping 


15.3x20.4 


1.3x1.8 


NGC 4388 


00086 
30572*^ 


Staring 
Staring 






NGC 4501 


30291*^ 


Staring 








03269 


Mapping 


15.3x20.4 


2.7x3.7 


TOLOLO 1238-364 


30291"^ 


Staring 








03269 


Mapping 


15.3x20.4 


3.9x5.2 


NGC 4968 


30291^ 


Staring 








03269 


Mapping 


15.3x20.4 


3.4x4.6 


M-3-34-64 


30323 


Staring 






NGC 5135 


00086 
30572*^ 


Staring 
Staring 






NGC 5194 


00159 


Mapping 


20.4x35.7 


0.90x1.6 


NGC 5347 


00086 
30572*^ 


Staring 

Staring 
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Table 5 — Continued 



Galaxy 


Program ID 


Mode 


Extraction Area^ 




57 


kpc 


Mrk 463 


00105 


Staring 






NGC 5506 


00086 
30572^ 


Staring 
Staring 






NGC 5929 


20140^ 


Staring 








03269 


Mapping 


15.3x20.4 


2.7x3.6 


NGC 5953 


00059 


Staring 








03269 


Mapping 


15.3x20.4 


2.2x2.9 


Arp 220 


00105 


Staring 






NGC 6890 


30291^ 


Staring 








03269 


Mapping 


15.3x20.4 


2.4x3.2 


IC 5063 


00086 
30572^ 


Staring 
Staring 






NGC 7130 


00086 
30572^ 


Staring 
Staring 






NGC 7172 


00086 
30572^ 


Staring 
Staring 






NGC 7582 


00024 
40936^ 


Staring 
Staring 






NGC 7590 


30291^ 


Staring 








03269 


Mapping 


15.3x20.4 


1.4x1.9 


NGC 7674 


03605 


Staring 








03269 


Mapping 


15.3x20.4 


8.5x11.4 



^The Program ID for the Sy2s in the [OIII] sample is 30773. 

^For low resolution spectral mapping observations. 

•^Galaxies that had dedicated off-source high-resolution ob- 
servations and were therefore background subtracted. See Sec- 
tion 3.2 for details. 
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Table 6: MIR Flux and PAH EW values for [OIII] Sample 



Galaxy 


[NeII]12.8lAim 


[OIV]25.89Aim 


Fmir 


EW 


^20-30^™ 










11.3/.jm 


17/im 






lO-i-* erg 


s~-^ cm~^ 


10~" erg s-^ cni-2 


fim 


/im 




NCG 0291 


20.6±1.2 


26.9±0.7 


4.87 


0.823 


0.471 


2.76 


Mrk 0609 


19.9±3.5 


7.82±1.20 


5.64 


1.021 


0.831 


2.22 


IC 0486 


4.20±0.83 


11.6±0.7 


4.59 


0.139 


0.147 


1.51 


2MASX J08035923+2345201 


0.96±0.71 


3.18±0.51 


0.63 


0.232 


0.552 


0.44 


2MASX J08244333+2959238 


3.07±1.08 


9.87±1.08 


8.36 


0.044 


0.094 


0.74 


CGCG 064-017 


2.52±0.27 


4.44±2.59 


0.73 


0.458 


0.373 


1.56 


2MASX J10181928+3722419 


0.71±0.69 


2.82±1.4 


1.10 


0.163 


0.107 


1.22 


2MASX J11110693+0228477 


1.15±0.68 


2.83±1.11 


1.56 


0.100 


0.272 


0.14 


CGCG 242-028 


1.59±0.77 


2.14±0.37 


0.39 


0.110 


0.675 


0.26 


SBS 1133-f572 


13.0±0.9 


16.1±3.8 


9.63 


0.135 


0.131 


1.57 


Mrk 1457 


9.87±0.89 


3.22±0.78 


2.93 


0.631 


0.421 


1.92 


2MASX J11570483-h5249036 


0.91±0.67 


5.53±0.93 


2.33 


0.054 


0.115 


1.52 


2MASX J12183945-h4706275 


<1.23 


4.31±0.18 


1.32 


0.048 


0.018 


0.76 


2MASX J12384342-h0927362 


<2.25 


6.46±1.13 


3.02 


0.007 


0.090 


0.34 


CGCG 218-007 


10.7±0.7 


17.0±0.8 


2.95 


0.887 


0.435 


2.66 


2MASX J13463217-h6423247 


2.10±1.19 


2.02±0.47 


0.86 


0.266 


0.183 


1.50 


NGC 5695 


2.94±0.58 


7.64±0.20 


0.87 


0.459 


0.666 


1.83 


2MASX J15522564+2753435 


- 


8.04±1.21 


1.44 


0.019 


0.059 


0.86 


SBS 1609-^490 


<1.53 


8.45±0.84 


2.35 


0.037 


0.098 


0.60 


2MASX J16164729-f3716209 


1.78±0.79 


5.11±0.54 


1.73 


0.051 


0.066 


1.38 
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Table 7: MIR Flux and PAH EW values for 12;mii Sample 



Galaxy [NeII]12.81/im [OIV]25.89/im Fmir EW 020-30^™ 

11.3/ini 17yum 

10~^^ erg s~^ cm~^ 10^^^ erg s^^ cm^^ /im /im 

18.2 0.058 0.130 -0.08 

4.28 0.347 1.92 

8.09 0.135 0.117 0.51 

6.59 0.131 0.239 1.57 

10.3 0.131 1.00 
2.76 0.578 1.35 1.45 
7.98 0.129 2.74 
4.19 1.85 4.13 
2.39 1.63 0.933 1.53 

6.44 0.083 0.187 1.79 
1.41 0.947 1.86 1.35 
11.7 0.253 0.065 0.82 

7.78 0.174 0.105 0.12 
11.7 0.010 0.035 0.85 

4.16 0.387 0.322 2.46 
8.07 1.71 1.63 2.06 

4.04 0.086 0.060 0.54 
6.72 0.016 0.58 

17.1 0.048 0.257 1.37 
0.58 1.15 0.792 1.61 

4.17 2.28 1.82 2.55 

7.45 5.512 6.27 

4.05 0.380 0.357 0.91 

14.2 0.009 0.063 0.97 

4.46 0.138 0.113 2.51 
3.13 0.170 2.31 2.53 

5.06 1.51 0.275 2.70 

1.79 1.95 1.95 1.80 
11.1 0.012 0.103 0.64 

^Program ID 30572 only for Fmir- 
^Program ID 00086 only for Fmir- 
^Program IDs 30572 and 00086 averaged together for Fmir- 



IRAS 00198-7926 


9 p;7+0.19 


25.7111 


NGC 0424 


q OQ+0.70 


OK 9+3.3 
^"^•^-3.0 


NGC 1068 


1491}? 


17631??! 


NGC 1144 


D-UU_o.5g 


4.911°:?^ 


NGC 1320 


9 CI +0.26 
^•^-■--0.20 


2i.6lS:° 


NGC 1386^ 


4 74+0.66 


82.01^:° 


F04385-0828 


A 90+0.42 


8.241}:^^ 


NGC 1667 


q qn+0.30 
'-'•'-"J-0.30 


2-O'-0.35 


F05189-2524 


r-89trel 


<69 


F08572+3915 


2.45l°:r7 


<33 


NGC 3982 


3.101°:}? 


4 r,9+0.52 


NGC 4388^ 


24.6111 


2621^ 


NGC 4501 


1.631°;}^ 


9 71+0.23 
^•' -'--0.21 


TOLOLO 1238-364 


13.91°:! 


14.9+1? 


NGC 4968 


8.081°:!^ 


26.3+f° 


M-3-34-64 


16 9+°-^° 
iD.y_0.75 


91.11-9^ 


NGC 51352 


32.6l|e^ 


56.71}?:^ 


NGC 5194 


22.9l°:« 


27.21?:^ 


NGC 53473 


1.601°:^^ 


6.391}:^^ 


Mrk 463 


2.99l°:i 


54.911? 


NGC 5506^ 


21.01}:^ 


2161}^ 


NGC 5929 


3.08+°-?^ 


4.381°:^! 


NGC 5953 


i7.8i°:« 


15.81}:? 


Arp 220 


21.81}:° 


<88 


NGC 6890 


3.481°:^° 


8.071°:^^ 


IC 5063 


7 0-1 +0.84 
'•Oi_o.8o 


97.0+?-^ 


NGC 7130^ 


22.31}:^ 


15.81?:° 


NGC 7172 


11.21}:^ 


34.9l}:e^ 


NGC 7582 


77.91^7 


87.lltl 


NGC 7590 


9 49+0.14 


3.611°:!^ 


NGC 7674 


R 04+0.49 


37.21^° 
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Table 8: Diagnostic Ratios: Optical & Infrared 



Diagnostic Ratio 


Sampl( 


n 




Mean 
dex 


dex 


^[OIII],corr/^[OIII],obs 


Both (51 Sy2s) 






0.73 


0.52 




[OIII]-only (20 Sy2s) 






0.34 


0.22 




12//m-only (31 Sy2s) 






0.97 


0.51 




Sy2s (Winter et al. 2010, 24) 


0.43 


0.49 




Syls (Winter et al. 2010, 29) 


0.28 


0.48 


'P[OIV]l'P[OIII],corr 


Both (48 Sy2s) 






-0.59 


0.48 




[OIII]-only (20 Sy2s) 






-0.28 


0.22 




12/im-only (28 Sy2s) 






-0.82 


0.49 




Sy2si (8) 






-0.38 


0.56 




Sylsi (10) 






-0.90 


0.58 


F[0/y]/F[0///],o6s 


Both (48 Sy2s) 






0.08 


0.41 




[OIII]-only (20 Sy2s) 






0.06 


0.29 




12/im-only (28 Sy2s) 






0.09 


0.48 




Sy2s (Melendez et al. 


2008, 


12) 


0.60 


0.74 




Sy2s (Diamond-Stanic 


et al. 


2009, 56) 


0.57 


0.67 




Syls (Melendez et al. 


2008, 


18) 


-0.21 


0.42 




Syls (Diamond-Stanic 


et al. 


2009, 16) 


-0.07 


0.53 


FAf//?/F[o///],corr 


Both (49 Sy2s) 






1.89 


0.63 




[OIII]-only (20 Sy2s) 






2.24 


0.26 




12/im-only (29 Sy2s) 






1.66 


0.70 


Fm7/?/F[0///],o6s 


Both (49 Sy2s) 






2.61 


0.63 




[OIII]-only (20 Sy2s) 






2.58 


0.37 




12//m-only (29 Sy2s) 






2.62 


0.77 




Syls2 (12) 






2.56 


0.50 


FM7/?/F[o/y] 


Both (46 Sy2s) 






2.46 


0.38 




[OIII]-only (20 Sy2s) 






2.52 


0.26 




12//m-only (26 Sy2s) 






2.41 


0.44 




Syls^ (24) 






2.59 


0.38 



Ratios in log space 

^Syl and Sy2 [OIII] values from Winter et al. (2010) and 
^Syl MIR values from Deo et al. (2009) and [OIII] values 
^Syl MIR values from Deo et al. (2009) and [OIV] values 
Melendez et al. (2008) 



[OIV] values from Weaver et al. (2010) 
from Melendez et al. (2008) 
from Tommasin et al. (2010) & 
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Table 9: Diagnostic Ratios: Optical & Infrared vs Radio^ & X-ray^ 



Diagnostic Ratio 



Sample 



Mean 
dex 



a 
dex 



F[OIII],corr/^8AGHz 
F[OIII],obs/FsAGHz 

F[OIV]/^SAGHz 

Fmir/^saghz 

Fl4- igSfceF /F[0///] ,1 



Fl4-195keV /F[OIII],obs 



^M-igbkeV /^[OIV] 



/F^ 



'U-ldbkeV/i" MIR 



Fl4-195fcey/F8.4Gi?2 



12/ini-only (26 Sy2s) 

12;Lini-only (26 Sy2s) 

Syls (12) 

12/xni-only (23 Sy2s) 

Syls (25) 

12/xm-only (25 Sy2s) 

Syls (22) 

Both (51 Sy2s) 

[OIII]-only (20 Sy2s) 

12/xm-only (31 Sy2s) 

Sy2s (Winter et al. 2010, 24) 

Syls (Winter et al. 2010, 29) 

Both (51 Sy2s) 

[OIII]-only (20 Sy2s) 

12/xni-only (31 Sy2s) 

Sy2s (Winter et al. 2010, 24) 

Syls (Winter et al. 2010, 29) 

Both (48 Sy2s) 

[OIII]-only (20 Sy2s) 

12/xm-only (28 Sy2s) 

Sy2s (Weaver et al. 2010, 33) 

Syls (Weaver et al. 2010, 37) 

Syls (Rigby et al. 2009, 17) 

Both (49 Sy2s) 

[OIII]-only (20 Sy2s) 

12;Lim-only (29 Sy2s) 

Syls3 (21) 

12/im-only (26 Sy2s) 

Syls (21) 



3.56 


0.87 


2.53 


0.99 


2.94 


0.93 


2.83 


0.64 


2.84 


0.90 


5.15 


0.66 


5.41 


0.78 


0.50 




1.79 




0.48 




2.30 


0.68 


1.74 


0.54 


1.48 




2.37 




1.53 




2.73 


0.57 


2.02 


0.54 


1.53 




2.14 




1.54 




2.30 


0.47 


2.60 


0.47 


2.68 


0.39 


-0.78 




-0.47 




-0.72 




-0.28 


0.26 


4.27 




5.30 


0.61 



Ratios in log space 

^Syl radio data from Thean et al. (2000) with [OIII] and [OIV] data from 

Melendez et al. (2008) and Diamond-Stanic et al. (2009), MIR values from 

Deo et al. (2009) and 14-195 keV data from Tueller et al. (2009) and 

Melendez et al. (2008). 

^Mean of Sy2 X-ray flux ratios calculated using survival 

analysis (Kaplan-Meier estimator) to include upper limits. 

The KM-estimator requires at least 2 detections, so one of the 

[OIII]-sample upper limits was converted to a detection which biases 

the mean for this sub-sample. 

^Syl MIR data from Deo et al (2009) with 14-195 keV values from 

Melendez et al. (2008) and Tueller et al. (2009). 
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Table 10: Results of Two Sample Tests between Syls and Sy2s 









Kuiper Test 




KS Test 


Diagnostic Ratio 


Neff 


D 


Probability 


D 


Probability 


F[o/y]/F[o///],ofes 


19.9 


0.310 


0.223 


0.301 


0.042 


FM/_R/F[o/77],obs 


8.49 


0.274 


0.912 


0.160 


0.949 


Fm/_r/F[o/t/] 


15.8 


0.264 


0.674 


0.203 


0.487 


F[OIII],obs/F8AGHz 


8.21 


0.385 


0.544 


0.340 


0.242 


F[OIV]/F8AGHz 


12.0 


0.318 


0.576 


0.200 


0.684 


Fmir/Fsaghz 


11.7 


0.344 


0.457 


0.253 


0.389 



Neff = ni X ^2/(^1 + 722) where ni and ^2 
are the number of data points in each sample. 



Table 11: Correlation of Diagnostic Ratios with Eddington Parameters 



Diagnostic Ratio 


P 


p 

-■- uncorr 


Optical & Infrared 


F[OIII],corr/F[oin],obs 


-0.224 


0.115 


F[OIV]/F[OIII],corr 


-0.131 


0.375 


F[OIV]/F[OIII],obs 


-0.351 


0.015 


FMIR/F[oin],corr 


-0.164 


0.259 


FMIR/F[oiII],obs 


-0.376 


0.008 


Fmir/F[oiv] 


-0.074 


0.625 


Radio and X-ray 


F[OIII],corr/FsAGHz 


0.008 


0.967 


F[OIII],obs/F8AGHz 


0.142 


0.491 


F[OIV]/F8AGHz 


-0.220 


0.314 


Fmir/Fsaghz 


-0.415 


0.039 


Fl4:-195keV /F[Oin],corr 




0.329 


Fu-ig5keV /F[oin],obs 




0.194 


Fi4-i95fcey/F[o/y] 




0.914 


Fl4-195fcey/FM/_R 




0.449 


Fl4~195fcey/F8. 40/^2 




0.283 



X-ray correlations computed using Cox Hazard Model 
in ASURV software package. 
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A. Aperture Bias 

As the optical data for tlie 12/im sample were culled from the literature, we examined 
the data to see if an aperture bias was evident. The aperture sizes used for the optical data 
(fioiii) ranged from 3 to 14" for the 12/im sources and was 3" for all the SDSS Sy2s. For 
several of our sources, we found the size of the NLR (^atlr) from Schmitt et al (2003a) and 
estimated the NLR for the remainder using log Rmaj = (0.31 ± 0.04) x log L[OIII] - 10.08 
± 1.80 (Schmitt et al 2003b). As shown in Figure IAT] (same color coding as Figure 1 in the 
main text), the aperture was large enough to encompass the full NLR for all sources. Hence, 
we are not "missing" any of the NLR optical flux. 

But are we observing too much [OIII] flux, perhaps from starburst contamination in the 
host galaxy which would affect [OIII] emission more than [01 V]? If this is the case, we would 
expect the F[oiv]/Fioiii],obs ratio to decrease as the projected size of the aperture increases 
and the "PAH-strong" sources to lie at systematically higher aperture sizes. However, neither 
of these trends are apparent (Figure rA.2p . reaffirming the results in the main text where we 
find no evidence for starburst bias on [OIII] emission. We also note that the opposite effect 
is also absent, namely increase of F[oiv]/Fioiii],obs with aperture size. This indicates that 
though the [OIV] flux is collected from the Long-High module which has a larger aperture 
than the [OIII] data, the [OIV] flux is not produced in regions outside of the NLR. Hence 
the dispersion present in the flux ratios is not due to sampling the galaxies at different scales 
between the separate IR and optical observations. 



B. Starburst Contribution to the MIR 

As discussed in the main text, in many cases the Spitzer IRS data show evidence for the 
presence of both an AGN and active star formation. This implies that the MIR continuum 
will include emission from warm dust heated by both the AGN and massive stars. Since we 
are interested in using the MIR luminosity as an indicator of the intrinsic luminosity of the 
AGN, we would need a way to subtract off the contribution from the dust heated by stars. 

We attempted to make this correction by using a simple dilution model based on the 
EW of the PAH features bracketing the continuum around 15.5/im (the PAH EW complexes 
at 11.3/im and 17/im, c.f. Genzel et al. 1998, Wu et al. 2009). This approach assumes that 
for starbursts there is a simple linear relationship between the PAH luminosity and that of 
the MIR continuum, and that the only effect of the AGN is to produce an additional source 
of MIR continuum emission, while not affecting the PAH luminosity. 

As we have shown, the PAH EWs are indeed correlated with other indicators of star 
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formation activity. We therefore compared the average PAH EWs for starburst galax- 
ies (< EWii,^^rn,SB >=2.87 and < EWn^^n^sB >=1.74, O'Dowd et al. 2009; Treyer 
private communication) to the measured EWs for our Sy2s and used this to derive the 
fraction of the MIR emission that is attributed to the AGN heating of the torus (i.e. 
Iagn = 1 — EWsy2/ < EWsB >)• In most cases, this correction was neghgible. How- 
ever, in the composite (strong PAH) AGN, this procedure seemed to subtract too much 
MIR flux, leading to poorer agreement with the other proxies and strong systematic trends 
with AGN luminosity (e.g. Figure R73] as compared with Figure [3]) • 

Such results indicate that accurately isolating the amount of MIR flux due to the AGN 
in strongly composite systems requires more detailed modeling that takes account of the 
possible influence of the AGN on the PAHs themselves. In the tight linear correlation we see 
between the uncorrected MIR and [OIV] luminosities (Figure 3) there is no evidence that 
the strong-PAH sources have systematically higher MIR luminosities. Taken at face value 
this would imply that in these sources the starburst does not dominate the MIR continuum. 
This would in turn imply that the EW of these long-wavelength PAHs with respect to the 
pure starburst MIR continuum is unusually high. This can not be understood as resulting 
from the destruction of PAHs by the AGN (e.g. Treyer et al. 2010), as the effect is in the 
opposite direction. 
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0.0 0.2 0.4 0.6 

^nlr/^oiii 



0.8 



1.0 



Fig. A.l. — Ratio of the NLR size to the optical aperture size used for the [OIII] flux and 
Ha and H/3 values. In all cases, the optical aperture encompasses the full NLR. 
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Fig. A. 2.— Left: Log {Fioiv]/F[oiii],obs) vs log {doiii)- Right: Log {F[oiv\/F[oiii\,obs) vs 
log {Onlr/Goiii)- The lack of significant trends suggest that the [OIV] flux originates in the 
NLR, consistent with the [OIII] flux, and that star formation processes in the host galaxy 
are not enhancing [OIII] emission. 
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Fig. A. 3. — Left: Distribution of F mir^agn lF{oiv\-, where Fmir,agn subtractions out the 
starburst contribution using the simple dilution model in the text. Right: Log {^ mir,agn) 
vs log (F[o/y]). Compared with the raw MIR flux, this corrected flux leads to a poorer 
agreement with the [OIV] flux and seems to over-substract the AGN contribution to the 
MIR flux for the "PAH-strong" sources. 



